Field Dependent Coherence Length in the Superclean, High-^ Superconductor 

CeCoIn 5 
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Using small-angle neutron scattering, we have studied the flux-line lattice (FLL) in superconduct- 
ing CeCoIns. The FLL is found to undergo a first-order symmetry and reorientation transition at 
~ 0.55 T at 50 mK. The FLL form factor in this material is found to be independent of the applied 
magnetic field, in striking contrast to the exponential decrease usually observed in superconductors. 
This result is consistent with a strongly field-dependent coherence length in CeCoIn 5 , in agreement 
with recent theoretical predictions for superclean, high-re superconductors. 

PACS numbers: 74.25.Qt, 74.25.Op, 74.70.Tx, 61.12.Ex 



Since the discovery of heavy-fermion superconductiv- 
ity in CeCoIns, a plethora of interesting phenomena have 
been observed in this material. Among these are one 
of the highest critical temperature (T c = 2.3 K) in any 
heavy-fermion superconductor Q, ci-wave pairing sym- 
metry field- and pressure- induced quantum-critical 
points and non- Fermi liquid behavior 0,0,0,0,01, strong 
indications of the first realization of a non- uniform Fulde- 
Ferrell-Larkin-Ovchinnikov state H E, E E E , 
and suggestions of multi-band or multi-order parameter 
superconductivity EjEI- 

In addition, CeCoIns is also found to represent an ex- 
treme case of a clean, high-re superconductor. The elastic 
electronic mean free path in this material ranges from / 
= 840 A at T c El , increasing exponentially as the tem- 
perature is decreased and reaching values of 4 to 5 fj,xn 
at 400 mK EE. Literature values for the penetration 
depth varies from 190 nm El t° 280 nm El- Estimates 
of the orbital critical field based on dH^/dT^,. range 
from i?J r c b (0) = 13.2 T0to 15.0 T E, yielding an in- 
plane coherence length of £ rb = 47 — 50 A, and hence re 
= 40 - 60 and l/£ ~ 10 3 or larger at temperatures below 
a few hundred milliKelvin. 

In many superconductors, detailed information about 
the nature of the superconducting state has been ob- 
tained from small-angle neutron scattering (SANS) stud- 
ies of the flux- line lattice (FLL). Examples of this are 
FLL symmetry transitions driven by non-local electrody- 
namics E or a superconducting gap anisotropy E>El' 
and subtle changes in the fundamental length scales ob- 
tained from the FLL neutron reflectivity E, • I n this 
Letter we report on SANS studies of the FLL in CeCoIn 5 . 
In striking contrast to the exponential dependence usu- 
ally observed in superconductors, the FLL form factor 
in this material is found to be independent of the ap- 
plied magnetic field. This result indicates a coherence 



length which decreases with increasing field, in quali- 
tative agreement with recent theoretical predictions for 
clean, high-re superconductors |27j |. 

The SANS experiments were carried out at the D22 
and Dll instruments at the Institut Laue-Langevin. The 
CeCoIn 5 single crystals used in the experiment were 
grown from excess indium flux, and had a T c = 2.3 K 
and a H C 2(0) = 5.0 T for fields applied parallel to the c 
axis . The sample was composed of three individually 
aligned single crystals with thicknesses 0.13 - 0.22 mm 
mounted side by side. The total mass of the sample was 
86 mg. The use of rather thin crystals was necessary, due 
to the strong neutron absorption by indium. Incident 
neutrons with a wavelength of A n = 4.5 A (Dll) and 
7 A (D22) and a wavelength spread of AA„/A n = 10% 
were used, and the FLL diffraction pattern was collected 
by a position sensitive detector. For all measurements, 
the sample was field-cooled to a base temperature of 40 
- 50 mK in a dilution refrigerator insert, placed in a su- 
perconducting cryomagnet. Horizontal magnetic fields in 
the range 0.4 - 2.0 T were applied parallel to the crys- 
talline c axis and the incoming neutrons. Background 
subtraction was performed using measurements follow- 
ing zero-field cooling. 

Fig. ^shows FLL diffraction patterns obtained at three 
different applied fields. Each image is a sum of the scat- 
tering from the FLL, as the sample is rotated and tilted 
in order to satisfy the Bragg condition for the different 
reflections. At fields below 0.55 T, 12 (2 x 6) reflections 
are observed as shown in Fig. ^a), corresponding to 
two nearly hexagonal domains with Bragg peaks aligned 
along (lOO)-directions. As the field is increased above 
~ 0.55 T the FLL undergoes a first order transition to a 
rhombic symmetry as shown in Figs. Q[b) and (c). Again 
two rhombic FLL domain orientations are observed, in- 
dicated by the 8 (2 x 4) Bragg peaks. As evident from 
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FIG. 1: (Color) FLL diffraction patterns for CeCoIn 5 with applied fields 0.5 T (a), 0.55 T (b) and 0.75 T (c), after subtraction 
of background measurement. The data is smoothed and the center of the image is masked off. The crystalline a axis is vertical. 



the decreasing peak splitting in the rhombic phase, the 
FLL gradually evolves towards a square symmetry as the 
field is increased. These results are in agreement with our 
earlier studies [28j . 

The evolution of the symmetry transition, quantified 
by the FLL opening angle /3, is summarized in Fig. [21 
Though it was not possible to reliably fit and extract a 
FLL split angle above 0.85 T, high-resolution measure- 
ments up to 1.0 T still showed a weak rhombic distortion. 
A linear extrapolation of the opening angle to (3 = 90°, 
yields a transition field H2 = 1.1 T above which a per- 
fectly square FLL is realized. Besides the opening angle, 
we also used scattering vectors belonging to the same do- 
main to obtain the internal field by the relation q\ x q-i = 
(2ir) 2 B/(/) 0j where (f> = 20.7 x 10 4 TA 2 is the flux quan- 
tum. From this we obtain dB/d(fi H) = 0.992 ± 0.012, 
which leads us to set B = /j, H for the remainder of this 
work. 

A square FLL can cither be stabilized by a gap 
anisotropy as observed e.g. in YBCO , or by non-local 
electrodynamics coupled with a Fermi surface anisotropy 
as seen in the borocarbide superconductors [22I l29j . In 
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FIG. 2: Field dependence of the FLL opening angle, (3. The 
insets show the real space orientation of the FLL unit cell in 
the low and high field configurations. 



the case of CcCoIn 5 the orientation of the gap nodes have 
been subject to controversy 0, 0, 0], although recent 
theoretical work aimed at resolving this issue concluded 
that the pairing symmetry in this material is d x 2_ y 2 |30l| . 
As we have previously reported, such an orientation of 
the gap nodes is consistent with the high-field square FLL 
being stabilized by e?-wavc pairing |2£(. On the other 
hand, an extrapolation to H = of the the high-field 
opening angle in Fig. ^yields (5 « 60°, as expected if the 
FLL symmetry is determined by non-local effects |29j . 

We now turn to measurements of the FLL form factor, 
which are the main results of this Letter. Fig|3{a) shows 
the FLL reflectivity obtained from the integrated inten- 
sity of the Bragg peaks, as the sample is rotated through 
the diffraction condition. The reflectivity is given by 



(1) 



where 7 = 1.91 is the neutron gyromagnetic ratio, A n is 
the neutron wavelength, t is the sample thickness, and q 
is the scattering vector j^lj. The form factor, F(q), is the 
Fourier transform of the magnetic field profile around a 
vortex, and depends on both the penetration depth and 
the coherence length. Fig. |3Jb) shows the FLL form fac- 
tor determined from the reflectivity. Here we have used 
the measured opening angle, @{H), in Fig. |2]to determine 
the magnitude of the scattering vector, and hence com- 
pensate for effects due to the FLL symmetry transition. 
In the high-field rhombic FLL phase q sq = q Q j \/sin /3 
where q Q = B / <p . In the low-field distorted hexag- 
onal phase, the 4 Bragg peaks aligned along the (100)- 
direction have ghexi = q \f2{l — cos/3)/ sin/3, while for 
the remaining 8 peaks ghex2 = qo / \/ sin /3. 

The field-independent form factor observed for 
CeCoIns is in striking contrast to the exponential de- 
crease observed in other superconductors. However, it is 
important to note that if the form factor was truly inde- 
pendent of q, this would imply an unphysical situation 
with a diverging magnetic field at the vortex center and 
a coherence length equal to zero. To reconcile this appar- 
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FIG. 3: Field dependence of the reflectivity (a) and form 
factor (b) for the (1,0) FLL Bragg reflections. Triangles cor- 
respond to a hexagonal FLL and the squares and diamonds 
to a rhombic FLL. The solid square is from previous work 
|2Sj . The dashed and dotted lines are calculated using Eqs. 
Q and with the following values for the coherence length 
and the penetration depth: £ = 81 A, A = 2350 A (dashed 
line); £ = 50 A, A = 3750 A (dotted line). The solid line 
corresponds to a constant form factor, F = 2.08 x 10 -4 T. 



ent paradox, it is necessary to assume a field dependence 
of either the superconducting penetration depth, the co- 
herence length or both. 

Several models exists for the magnetic field profile 
around the vortices. In the following we have based the 
analysis on the form factor obtained by Clem 



F(q) = B 



9l<i(g) 
1 + AV 
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V2e(q 2 +X~ 2 ) 1/2 , (2) 



valid for k 1. It is worth to point out that while this 
form factor was obtained for an isotropic superconduc- 
tor, numerical simulations found no significant difference 
between superconductors with s- or d-wave pairing [33j |. 

The dashed lines in Figs. Efa) and (b) show the re- 
flectivity and form factor calculated from Eqs. and 
@ using an average literature value for the penetration 

depth A = 2350 A and £ c2 = J ^f-^- = 81 A. As al- 
ready discussed, these calculated values are not consis- 
tent with the experimental data. However, a more ap- 
propriate value for the coherence length is obtained from 

the orbital critical field, £ orb = J 2 j£ « 50 A [l7l l2l| . 



This does correspond to a smaller slope of |F((7)| 2 , but 
also increases the discrepancy between the measured and 
calculated values of |-F(<?)| 2 since, for a given £ and q (or 
H), the magnitude of the form factor is determined by 
A 2 . Using the penetration depth as a fitting parameter, 
one obtains the dotted lines in Fig. [21 for A = 3750 A. 
Although this provides a better agreement with the data, 
such a large value of A is not consistent with reports in 
the literature 0,H3|- I n contrast, a perfect fit to the 
data is obtained with a constant value of the form factor, 
F = 2.08 x 10~ 4 T, and correspondingly a reflectivity 
cx 1/q, as shown by the solid lines in Fig. [3J 

It is important to note that no significant disordering 
of the FLL is observed. Except for systematic differ- 
ences between the two beamlines, the FLL rocking curve 
widths remain essentially constant throughout the mea- 
sured field range. On D22 we find rocking curve widths 
going from 0.14°±0.02° FWHM at 0.55 T to 0.16°±0.01° 
FWHM at 1 T, comparable to the experimental resolu- 
tion (~ 0.08° FWHM). Below the reorientation transi- 
tion a slightly higher value of 0.20° ± 0.04° FWHM is ob- 
served. On Dll FLL rocking curve widths decrease from 
0.26° ± 0.03° FWHM at 0.7 T to 0.19° ± 0.01° FWHM 
at 2 T. Such narrow rocking curve widths indicate a very 
well ordered FLL with a longitudinal correlation length 
in the micron range, consistent with weak pinning due to 
the high cleanliness of CeCoIns. Furthermore, FLL dis- 
order above a certain threshold has been shown to lead 
to decrease in the scattered intensity, exceeding the usual 
exponential field dependence of the form factor j^j • We 
can therefore exclude FLL disordering as the cause for 
the field- independence of the form factor in CeCoIns. 

A constant form factor could in principle be due to a 
decrease of the penetration depth with increasing field, 
caused either by an increasing superfluid density or a 
non-uniform spin magnetization contributing to the mag- 
netic flux carried by each vortex. We do not consider the 
first possibility realistic. Furthermore, while CeCoIns is 
indeed paramagnetic [2lj|, spin polarization effects will 
lead to an enhancement of the form factor in contrast to 
the strong suppression observed in this material [35t l3(| . 
We therefore conclude that while paramagnetic effects 
may contribute, they are not the dominating mechanism 
behind the field-independent form factor. In the follow- 
ing we therefore restrict our analysis to consider only a 
field-dependent coherence length. 

In Fig. QJ we show the coherence length obtained by 
varying £ to achieve the measured form factor at each 
field. The coherence length is found to follow a 1/yH 
behavior. While different models for F(q) will provide 
slightly different values for the coherence length, the 
qualitative behavior will remain unchanged. Fig. [3] also 
shows the extracted coherence length as a function of 
intervortex spacing, a. Within the experimental error 
shown by the scatter in the data, the coherence length 
is found to increase linearly with the intervortex spacing 
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FIG. 4: Field dependence of the coherence length, (,(H) plot- 
ted versus applied field (solid symbols) , and vortex separation 
(open symbols). The solid and dashed lines are fits to the data 
described in the text. 



while at all times satisfying £ < a. Extrapolating to a 
vortex separation of 125 A (corresponding to an orbital 
critical field H orh = 13.2 T) yields £ = 43.0 ± 4.6 A, in 
reasonable agreement with £ or b = 50 A. 

We believe that our data provide the first example 
of a strongly field-dependent coherence length in a su- 
perclean, high-ft superconductor as recently theoretically 
predicted by Kogan and Zhelezina j^- In their model, 
the coherence length is predicted to be proportional to 
l/y/H and correspondingly depends linearly on the vor- 
tex separation. Experimentally we find d^/da = (0.55 ± 
0.02)/V27r, which is roughly twice the theoretically pre- 
dicted slope [27]] • However, the theoretical prediction was 
based on assumptions of weak-coupling s-wavc supercon- 
ductivity and a simple Fermi surface topology (sphere or 
cylinder), which are not valid for CeCoIn 5 . Based on this 
we argue that the agreement between theory and exper- 
iment is remarkable, supporting the universal nature of 
the effect. At low fields the model predicts the coherence 
length to saturate. Based on our experimental results 
this cut-off in CeCoIns occur below the lowest measured 
field of 0.4 T. We believe that the field-dependent coher- 
ence length is so prominently observed in CeCoIns due to 
the combination of a large k and the very high cleanliness 
of this material. 

Additional evidence for the unusual magnetic-field re- 
sponse of CeCoIns has been observed in the quasiparticlc 
mean free pat h, as extracted from measurements of the 
Hall angle |3j|. Specifically, the mean free path is found 
to decrease as the applied field increases, being roughly 
equal to the vortex separation for the range of fields cov- 
ered by our SANS measurements. This leads us to spec- 
ulate that a connection exists between £ and / beyond 
the simple Pippard model. 

Finally, we want to note that while previously reports 
of a field dependent vortex core size have been made by 
Sonier et al. 133, 1391 , these were attributed to a vortex- 



squeezing effect. However, since the variation of the vor- 
tex core size in these cases was significantly smaller than 
for CeCoIn 5 , we do not believe that such effects can ex- 
plain the data reported here. 

In summary, we have presented SANS measurements 
of the FLL in CeCoIns. These measurements indicate a 
strong reduction of the superconducting coherence length 
with increasing field, in qualitative agreement with recent 
theoretical predictions for superclean, high-K supercon- 
ductors. 
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